Contrasting population genetic structure of three semi‐terrestrial brachyuran crabs on the coast of the Japanese archipelago

Abstract Anthropogenic activities have reduced ecotones between the ocean and land, which is likely to threaten the population of brackish‐water brachyuran crabs. To assess the current status of these crabs, we examine the population genetic structures of three semi‐terrestrial brachyuran crabs widely distributed along the coast of the Japan and to clarify factors determining their genetic structures. We collected 184 Orisarma dehaani, 252 Chiromantes haematocheir, and 151 Helice tridens crabs from 36 localities of the Japanese archipelago. Genome‐wide SNP data from these crabs were analyzed using MIG‐seq. Bayesian clustering of STRUCTURE and DAPC analysis were used to identify genetically disturbed populations and to visualize genetic differentiation between local populations. Genetic population structure showed clear differentiation between populations on the Pacific coast of the Tohoku region and on other Japanese coasts in O. dehaani, but not in C. haematocheir or H. tridens. The inbreeding coefficient of O. dehaani was significantly higher on the Pacific coast of the Tohoku region compared to other Japanese coasts. C. haematocheir and H. tridens had homogeneous genetic structures along the Japanese coast, but showed genetic differentiation of a local population at their range limits. Thus, O. dehaani showed little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to ocean currents. Although such a regional differentiation was not found in C. haematocheir and H. tridens, one population of C. haematocheir was genetically isolated at the edge of its distribution range and likely vulnerable to environmental changes.


| INTRODUC TI ON
In recent decades, the development of coastal regions by human activities has altered and reduced the ecotone areas between the ocean and land.These areas are the habitats of brackish-water and semi-terrestrial brachyuran crabs (Miyake, 1983).Accordingly, the populations of these crabs may be threatened with extinction.Therefore, it is necessary to assess the current status of the brackish-water brachyuran crab populations.
Studies have shown that endangered populations tend to have low genetic diversity and become genetically isolated from populations in other areas (Frankham et al., 2010).Indeed, in addition to anthropogenic activities, various factors, such as the dispersal ability of the planktonic larval stages (DeWoody & Avise, 2000;Palumbi, 1992) and ocean currents as invisible physical barriers (Gilg & Hilbish, 2003;Palumbi, 1992), are known to affect the population genetic structure of coastal invertebrates.In the Japanese archipelago, several ocean currents, including Oyashio, Kuroshio and Tsushima currents, are known to genetically differentiate the local populations of coastal organisms, such as gobies (Matsui, 2022), snails (Kojima et al., 2004;Nakano et al., 2010;Yamakawa & Imai, 2014) and the seaweeds (Hu et al., 2016;Zhong et al., 2020).These organisms have phylogeographic patterns composed of "Pacific Ocean lineage" and "Sea of Japan lineage," which correspond to the paths of the two currents, the Kuroshio Current and the Tsushima Current, respectively (Matsui, 2022).Although several studies have examined the genetic population structure of brachyuran crabs, such as the varunid crab Ptychognathus ishii (Kawane & Wada, 2015), the sesarmid crab Clistocoeloma sinense (Yuhara et al., 2014), the camptandriid crabs Deiratonotus cristatus (Kawamoto et al., 2012) and Deiratonotus kaorie (Kawane et al., 2012), and the ocypodid crabs Austruca lactea (Tokuyama et al., 2020) and Tubuca arcuata (Aoki et al., 2008), the factors affecting their genetic structure are not well understood, since the distribution ranges of these crabs are geographically restricted to the limited area of the Japanese coast.
The sesarmid crabs Orisarma dehaani and Chiromantes haematocheir, and the varunid crab Helice tridens inhabit the coastal land-ocean ecotone and are endemic to East Asia (Japan, Korea, China, and Taiwan) (Miyake, 1983;Sakai, 1976;Sakai et al., 2006;Shih, 2007).These semi-terrestrial crabs have planktonic zoeal larval stages (Cuesta et al., 2006) that can spatially disperse using ocean currents to expand their range in coastal areas.Adult individuals of these three crab species live primarily in estuaries and salt marshes.
However, their habitats differ slightly each other (Miyake, 1983;Sakai, 1976): C. haematocheir inhabits mainly inland forest areas, O. dehaani occurs mainly in areas close to water, such as wetlands, and H. tridens inhabits muddy shores of the mid to high intertidal zone.In addition, these species are known to have different settlement patterns, presumably due to different preference for salinity and temperature (Matsumoto et al., 2020;Saigusa, 1981).
In East Asia, the Japanese archipelago is surrounded by coasts on all sides and has a very long coastline of about 30,000 km (Fujikura et al., 2010).As these semi-terrestrial crabs are widely distributed along the coastline of the Japanese archipelago (Miyake, 1983), they are ideal organisms to study the factors affecting the population genetic structure of brachyuran crabs because these crabs, which inhabit in the land-ocean ecotone, are often classified as locally endangered species in Japan (Japanese Association of Benthology, 2012).In this study, therefore, we examined the population genetic structures of these three crabs in the Japanese archipelago to clarify whether they are genetically differentiated among the geographically distinct areas and, if so, what factors caused these genetic differentiations.Specifically, we analyzed samples collected from 36 sites in the Japanese archipelago using sequence data of mtDNA COI sequence and multiplexed inter-simple sequence repeat (ISSR) SNP data.We then used these data to assess whether there were local populations threatened with extinction.

| Specimens
The three semi-terrestrial crabs were collected as whole bodies by hand-digging in 2019 and 2020 from a total of 36 tidal flats and salt marshes across the four major islands of the Japanese archipelago (Figure S1, Table S1).We collected a total of 184 individuals of the sesarmid crab Orisarma dehaani from 21 of 43 localities and 252 individuals of the sesarmid crab Chiromantes haematocheir from 26 localities.We also collected a total of 151 individuals of the varunid crab Helice tridens from 16 localities.These specimens were either frozen at −20°C or stored in 95% ethanol until DNA extraction.

| DNA extraction
Genomic DNA was extracted using a DNeasy Blood & Tissue Kit (QIAGEN).The concentration of gDNA was measured using a "NanoDrop™ Lite" spectrophotometer (Thermo Scientific, Waltham, MA) and was adjusted to 1 ng/mL for PCR amplification.
All amplified samples were purified using ExoSAP-IT® (Affymetrix) and subsequent product sequencing was performed at Eurofins Genomics.The obtained DNA sequences were deposited in DDBJ (accession numbers LC768188 -LC768321).

| MIG-seq
Genome-wide SNP data were analyzed using the MIG-seq method (Suyama & Matsuki, 2015) using DNA samples from 1 to 10 individuals collected at each location in each species, i.e., 154 individuals collected at 18 locations for O. dehaani, 183 individuals collected at 22 locations for C. haematocheir, and 145 individuals collected at 16 locations for H. tridens.A MIG-seq library was constructed by two step PCR as described by Suyama et al. (2022).
Inter-simple sequence repeat (ISSR) regions were amplified by first PCR using MIG-seq primer set-1 (Suyama & Matsuki, 2015).The PCR products were checked using a Microchip Electrophoresis System (MultiNA, Shimadzu, Japan) with a DNA2500 Reagent Kit (Shimadzu), which was also used to quantify the second-round PCR products.In the subsequent steps for a second PCR, the complementary sequences for the binding sites of Illumina sequencing flow cell and the indices (barcodes) for each sample were added to the first round PCR products.The second PCR products were mixed and extracted ≥350 bp fragments.Sequencing was performed by Illumina MiSeq sequencers (MiSeq Control Software ver.2.0.12;Illumina, San Diego, CA) using a MiSeq Reagent Kit ver.Extremely short reads containing adapter sequences and low-quality reads were removed using Trimmomatic 0.39 (Bolger et al., 2014)  between stacks (M) = 2, maximum mismatches between loci when building the catalog (n) = 2. SNP selection was performed using populations program of Stacks.We selected only SNPs with more than 60% coverage among samples from each crab for MIG-seq.
In addition, SNPs where one of two alleles had less than 5% was filtered out.SNPs with high heterozygosity (Ho ≥ 0.6) were also removed.To avoid linked SNP, we used only first SNP per locus.

| Data analysis
The number of mtDNA haplotypes, and the haplotype-and nucleotide-diversity were calculated for each population of the three semi-terrestrial crabs using the computer software programs Arlequin ver.2.001 (Schneider et al., 2000) and DnaSP 5.10.1 (Librado & Rozas, 2009).Phylogenetic relationships among the haplotypes were inferred for each species using the statistical parsimony software TCS Networks (Clement et al., 2002) implemented in the package PopART (http:// popart.otago.ac.nz).The alignment data were made with the default settings in MUSCLE implemented in MEGAX (Kumar et al., 2018).
To identify genetically distinct populations, we performed Bayesian clustering using STRUCTURE ver.2.3.4 software (Pritchard et al., 2000), which assigns individuals into K clusters.The population structure was simulated with values of K = 1-10 under an admixture model, i.e., the correlated allele frequency model (Hubisz et al., 2009).All runs included 100,000 Markov chain Monte Carlo generations after a burn-in period of 100,000 iterations.Ten runs were performed for each value of K.The optimal K value was determined using the ΔK methods of Evanno et al. (2005) in STRUCTURE Harvester program (Earl & Vonholdt, 2012, http:// taylo r0.biolo gy. ucla.edu/ struc tureH arves ter/ ).Detection of genetically distinct populations was done by the analysis of SNP data using STRUCTURE software (Pritchard et al., 2000).For species that showed clear genetic differentiation between populations (optimal K > 1) in the STRUCTURE, we performed a discriminant analysis of principal components (DAPC; Jombart et al., 2010) to visualize betweenpopulation differentiation in the SNPs.This analysis extracts information from the data by first performing a principal component analysis on user-defined populations.Then, using the PCA factors as variables, a discriminant analysis is performed to maximize the inter-population component of variation.In the present study, the DAPC was implemented using the optimized number of principal components as determined by the a-score and conducted using the "adegenet" package (Jombart et al., 2010) for the statistical platform R 4.03 (R Core Team, 2020) after the text data file was converted to the Genepop file format.The number of polymorphic SNPs observed, and the expected heterozygosity for each population were calculated using POPULATIONS in Stacks (Catchen et al., 2013;Rochette et al., 2019).Fixation index (F IS ) for each population was calculated and their significance was tested by 1000 permutations using FSTAT 2.9.4 (Goudet, 2003).Phylogenetic relationships among the populations of O. dehaani were assessed by constructing a neighbor-joining network based on DA distances between the populations (Nei et al., 1983) using SplitTree 4 (Bryant & Moulton, 2004).
Statistically significant difference nodes were evaluated from bootstrap probabilities based on 1000 replicates of the neighbor-joining method using POPTREE2 (Takezaki et al., 2010).

| Estimation of spatial genetic structure and migration patterns
Using the SNP data obtained for each of the three semi-terrestrial crabs, pairwise population F ST values were estimated by the AMOVA using the GenAlEx 6.5 (Peakall & Smouse, 2012).Since missing data can be particularly problematic for pairwise distance-based analyses, large numbers of missing data should be minimized (Blyton & Flanagan, 2012).For this analysis, we used a dataset consisting of 36 loci in O. dehaani, 19 loci in C. haematocheir, and 47 loci in H. tridens with 90%, 90% and 95% sharing sites, respectively.Statistical significance levels for all pairwise tests were set at 0.05 after adjusting for multiple comparisons using the Bonferroni correction (Rice, 1989) in this analysis.To evaluate genetic differentiation at the population level, we also measured the geographic coastal distance between each sampling locality along the coastline using distances between Japanese seaports (The Japan Shipping Exchange, Inc., 1996).For each of three semi-terrestrial crabs, a Mantel test (Mantel, 1967) was performed to assess correlations between genetic differentiation (F ST /(1− F ST )) and geographic distance with 9999 permutations in the "ade4" package (Dray & Dufour, 2007) for R 4.03 (R Core Team, 2020).S2).Haplotype diversity values (h) and nucleotide diversity values (π) ranged from 0.607 to 1.000 and from 0.0022 to 0.0079, respectively.According to a statistically parsimonious haplotype network (Figure S2a), the dominant haplotype was Od1, and the other haplotypes were separated from Od1 by 1 to 7 mutational steps, and formed a "bush-like" shape (Figure S2a,b).1).The STRUCTURE analysis of the MIG-seq data, which supported the partitioning of O. dehaani into two genomic components (maximum ΔK = 201.570:at K = 2; Figure S3a), showed a clear distinction between 5 local populations on the Pacific coast of the Tohoku region and 13 local populations on other Japanese coasts (Figure 1a).The same result was obtained by the DAPC analysis, which showed that the first axis separated the Tohoku Pacific coast populations from the other populations (Figure 1b).In this analysis, the populations of these two regions separated from each other with 96% probability, indicating that O. dehaani were grouped into two phylogenetic groups (Figure S4).The relative directional migration rates from the Pacific coast of the Tohoku region (cluster 1) to other Japanese coasts (cluster 2), detected by STRUCTURE analysis, were 0.5, while the rate in the opposite direction was 1.0 (Figure S5).

| MIG-seq
Although the bootstrap analysis showed that the rate of the former diction was not significant (p > .05), the rate of the latter direction was significant (p < .001).

| Mitochondrial DNA
In this species, DNA segments of 499 bp in the mitochondrial COI gene were commonly sequenced with no insertions and no deletions.
In these segments, we found a total of 39 variable sites and identified 50 haplotypes.Twelve of the 50 haplotypes were shared among 2 to 19 localities, while the remaining 38 haplotypes were localityspecific (Table S3).Haplotype diversity values (h) and nucleotide diversity values (π) ranged from 0.533 to 0.917 and from 0.0015 to 0.0067, respectively.According to a statistically parsimonious haplotype network (Figure S6a), the dominant haplotype was Ch1, and the other haplotypes were separated from Ch1 by 1 to 6 mutational steps, and formed three "star-shaped" clusters (Figure S6a,b).at K = 2; Figure 2a; Figure S3b).The same result was obtained by the DAPC analysis, which showed overlap among the populations in the coasts of the Japanese archipelago (Figure 2b).

| Mitochondrial DNA
In this species, DNA segments of 490 bp in the mitochondrial COI gene were commonly sequenced with no insertions and no deletions.In these segments, we found a total of 30 variable sites different and 35 mtDNA haplotypes of H. tridens.Among them, 14 haplotypes were found more than 2 localities, while the other 21 haplotypes were locality-specific (Table S4).The haplotype diversity values (h) and nucleotide diversity values (π) ranged from 0.533 to 0.917 and from 0.0012 to 0.0047, respectively.According to a statistically parsimonious haplotype network (Figure S7a), the dominant haplotype was Ht1, and the other haplotypes were separated from Ht1 by 1 to 2 mutational steps, and formed three "star-shaped" clusters (Figure S7a,b).(range: 0.117 to 0.26; mean = 0.188) (Table 3).In the STRUCTURE analysis, the population did not genetically separate the populations (maximum ΔK = 38.365: at K = 2; Figure 3a; Figure S3c).This trend was confirmed by the DAPC analysis, which showed overlap between the populations (Figure 3b).

| A comparison of the spatial genetic structure of three semi-terrestrial crabs
Pairwise F ST comparisons of O. dehaani showed significant difference in the population genetics between the Pacific coast of the Tohoku region and other Japanese coasts (Table S5).For C. haematocheir, pairwise F ST values were significant between the southern Kyushu Island population (locality code KMNE) and the other populations (Table S6).The KMNE population showed significant F ST values against 7 of the 21 populations examined.The pairwise F ST values of H. tridens were not significant among the coastal populations of Japan, except for five pairs of populations (Table S7).
Pairwise genetic distance (F ST /(1 − F ST )) was significantly and positively correlated with the pairwise geographic distance in O.
We also found that the inbreeding coefficient (F IS ) was significantly higher in O. dehaani (0.326) than in C. haematocheir (0.238) and H. tridens (0.188) (Figure S9).Furthermore, the F IS of O. dehaani populations was significantly higher on the Pacific coast of the Tohoku region (F IS = 0.347) compared to the other Japanese coasts (F IS = 0.318) (Figure S10).

TA B L E 3
Genetic diversity parameters of the 16 H.tridens populations in the Japanese archipelago estimated from 1025 single nucleotide polymorphisms (SNPs).The present study showed that O. dehaani populations were genetically divided into two groups, i.e. those distributed along the Pacific coast of the Tohoku region (from FRK to ISM populations) and those distributed along other Japanese coasts.However, C. haematocheir and H. tridens populations did not show such a clear regional differentiation in genetic structure.Below, we first discuss why the spatial pattern of population genetic structure differed among these three species.We then discuss whether or not there are populations at risk of extinction.

| Effects of life history traits
The present study showed that, among the three semi-territorial crabs, O. dehaani had a significantly higher inbreeding coefficient than C. haematocheir and H. tridens.Many benthic marine organisms have pelagic larval stages that extend their range by drifting in ocean currents as plankters (Cowen & Sponaugle, 2009).In these organisms, the distance of the dispersal range is generally related to the duration of the planktonic larval stages (Levin, 2006;Shanks, 2009).
In C. haematocheir and H. tridens, zoeal stages I to V are planktonic, whereas in O. dehaani only zoeal stages I to IV are planktonic (Cuesta et al., 2006;Mia & Shokita, 1997;Terada, 1976).The number of planktonic zoeal stages seems to correspond to the pelagic larval period.In fact, Terada (1976) reported that the pelagic larval period need to migrate offshore and, therefore, remain in nearshore waters (Baba, 1978;Saigusa, 1981).In support of this, O. dehaani larvae are more tolerant to low salinities than C. haematocheir and H. tridens larvae, supporting the idea that they complete their life cycle without migrating offshore (Irawan & Kijima, 1994;Matsumoto et al., 2020).
If these reproductive and developmental traits limit the migratory range of the larvae, that may have driven the observation that O. dehaani populations are genetically differentiated into two biogeographically distinct groups.
However, the low dispersal ability of O. dehaani larvae cannot explain why the genetic structure was different between populations on the Pacific coast of the Tohoku region and those on other Japanese coasts.If their population genetic structure is determined by low dispersal ability alone, it is likely that the genetic structure is fragmented at different locations along the coast of the Japanese archipelago.Thus, their population genetic structure should be additionally shaped by other factors, such as ocean currents, which may limit larval transport and dispersal (Williams & Hastings, 2013).

| Effects of ocean current and temperature
It is well known that the Pacific Ocean from the Tohoku region to the east coast of the Boso Peninsula has a lower surface water temperature than other areas of Japan, except for the coasts of Hokkaido, due to the influence of the Oyashio Current (Nishimura, 1992).
Orisarma dehaani was not distributed along the coast of Hokkaido, which is strongly influenced by the Oyashio Current (Miyake, 1983).
Lake Jusanko on the Sea of Japan side in northwestern Tohoku is the northern distribution limit of this species (Yuhara & Suzuki, 2021; this study).On the Pacific coast side, we could not collect O. dehaani individuals at the sampling localities located at latitudes higher than Hirase & Ikeda, 2015).The Kuroshio Current acts as a conveyor for larvae of various tropical and subtropical brachyuran crabs to the high latitude regions of Japan (Nishimura, 1992).However, it leaves the coast at the Boso Peninsula and flows eastward as the Kuroshio Extension in central Japan at about 35° N (Qiu, 2001).Thus, this geographic pattern of the Kuroshio Current appears to functionally act as a barrier to genetically separate O. dehaani populations on the Pacific coast of Tohoku from those in other regions.Prefecture (Osawa, 2022;Schubart & Ng, 2020), and the Kumanoe River population (KMNE) is near the edge of its range.In general, peripheral populations exhibit low genetic diversity and greater genetic differentiation due to smaller effective population size and greater geographic isolation compared to geographically central populations (Eckert et al., 2008;Pandey & Rajora, 2012).In addition, genetic drift distinguishes the genetic structure of small populations from that of larger populations (Austerlitz et al., 1997;Signorile et al., 2014).Indeed, in Miyazaki Prefecture, where the KMNE population is located, C. haematocheir is listed as a "near-threatened species" due to its small population size (Miyazaki Prefecture, 2020).

| Local isolation
The present study confirmed such local caution, as C. haematocheir populations in this marginal area of their range are somewhat genetically isolated and likely vulnerable to environmental changes.Note that these semi-terrestrial crabs are also distributed in coastal areas of continental East Asia, which are located to the west and south of the Japanese archipelago.Therefore, it is possible that some of the local populations in the Japanese archipelago have the same genetic structure as the East Asian populations.However, the population genetic information of semi-terrestrial crabs on the continental coasts is still limited.In future studies, it is necessary to clarify how the semi-terrestrial crab populations in the Japanese Archipelago are genetically related to those in East Asia.

| CON CLUS ION
The present study revealed that O. dehaani shows little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to the influence of the Kuroshio and Oyashio currents.Thus, these Tohoku Pacific populations appear to be isolated and worthy of conservation as independently evolving population groups (i.e., Evolutionarily Significant Units, ESUs) (Moritz, 1999).In contrast, the genetic struc-

ACK N OWLED G EM ENTS
This study was financially supported by KAKENHI (19K20488 and 23H02548) and the Grant for Environmental Research Project (193046) by the Sumitomo Foundation.We also thank the following persons who provided or were helping us collect the important samples of three semi-terrestrial brachyuran crabs for this study: Takao Suzuki, Tomomichi Sudo, Masatoshi Matsumasa, Hirokazu Abe, Kotaro Kan, Kenji Hirashima, Kyoko Kinoshita, Takashi Uchino, Kazuyo Manabe, 3 (150 cycle; Illumina) with 80-bp paired-end reads.The MIG-seq data was deposited at the DDBJ Sequence Read Archive (DRA, accession number DRA016447).
Finally, the gene flow between the populations of O. dehaani was estimated by using divMigrate-Online (https:// popgen.shiny apps.io/ divMi grate -online/ ).This program generates a migration network graph with relative values for gene flow among populations scaled to the largest magnitude estimated.For the analysis, we selected two populations obtained from the K clusters of the STRUCTURE analysis.We used Nm as a measure of genetic distance.The significance of asymmetric gene flow between populations was tested by 1000 bootstrap iterations.3 | RE SULTS3.1 | Orisarma dehaani3.1.1 | Mitochondrial DNADNA segments of 498 bp in the mitochondrial COI gene were sequenced and used for analyses.No insertions or deletions were found in these segments.In Orisarma dehaani, 48 different mtDNA haplotypes were identified from a total of 50 variable sites in the segments.Sixteen of the 48 haplotypes were shared among 2 to 19 localities, while the remaining 29 haplotypes were specific to single localities (Table

F
I G U R E 1 Results of Bayesian analysis of population structure (STRUCTURE) with K = 2 clusters and discriminant analysis of principal components (DAPC) for 18 populations of O. dehaani.(a) STRUCTURE composition plot for each individual grouped by the sampling localities.The locality number and code from Table 1 are shown on the panels.(b) Plots of the first two axes obtained by DA.Individuals from each population are represented by a population-specific symbol.The number of PCs retained, as determined by the a-score, was 58.Eigenvalues of the PCA and DA are inserted.

F
I G U R E 2 Results of Bayesian analysis of population structure (STRUCTURE) with K = 2 clusters and discriminant analysis of principal components (DAPC) for 21 populations of C. haematocheir.(a) STRUCTURE composition plot for each individual grouped by the sampling localities.The locality number and code from Table 2 are shown on the panels.(b) Plots of the first two axes obtained by DA.Individuals from each population are represented by a population-specific symbol.The number of PCs retained, as determined by the a-score, was 41.Eigenvalues of the PCA and DA are inserted.
was 17 and 15 days for C. haematocheir and H. tridens, respectively, while it was 12 days for O. dehaani.This observation suggests that O. dehaani has the lowest dispersal ability among the three semiterrestrial crabs.Saigusa (1981) showed that C. haematocheir and H. tridens released larvae in a clear semi-annual rhythm, whereas the larval release activity of O. dehaani did not have such a clear semi-annual rhythm.The observations suggest that O. dehaani larvae do not

F
I G U R E 3 Results of Bayesian analysis of population structure (STRUCTURE) with K = 2 clusters and discriminant analysis of principal components (DAPC) for 16 populations of H. tridens.(a) STRUCTURE composition plot for each individual grouped by the sampling localities.The locality number and code from Table 3 are shown on the panels.(b) Plots of the first two axes obtained by DA.Individuals from each population are represented by a population-specific symbol.The number of PCs retained, as determined by the a-score, was 72.Eigenvalues of the PCA and DA are inserted.

Furukawa-
Numa Lagoon, i.e., NUS, ASZ, TKH, TKS, MYK, ORK, and OTM(Yuhara & Suzuki, 2021; this study).The result suggests that the northern limit of their range was limited by low temperature.Since the Pacific coasts of the Tohoku region are still under the influence of the Oyashio Current, O. dehaani populations in the Tohoku region may have adapted to relatively low temperature, as in the case of the goby fish Chaenogobius annularis and the abalone Haliotis discus, which have different genetic structures in the northeastern areas of the Japanese Pacific coast due to the influence of low seawater temperature originating from the Oyashio Current(Hirase et al., 2021;

For
C. haematocheir and H. tridens, neither mtDNA nor MIG-Seq analysis detected genetic separation between populations living along the coasts of the Japanese archipelago.However, C. haematocheir was genetically distinct from other local populations in the Kumanoe River population (KMNE) in Miyazaki Prefecture.One possible reason for this local genetic differentiation is the effect of population size at the edge of the range.The southern limit of the distribution of C. haematocheir in the Japanese archipelago is in Kagoshima tures of C. haematocheir and H. tridens were not affected by ocean currents and were spatially homogeneous along the entire Japanese coast.However, for C. haematocheir, genetic differentiation was observed in a local population at the southern limit of its distribution range.Since its population is vulnerable due to its genetically small population size, it is necessary to continue monitoring changes in its population size in order to conserve a local population.AUTH O R CONTR I B UTI O N S Takeshi Yuhara: Conceptualization (lead); formal analysis (equal); funding acquisition (lead); investigation (equal); methodology (equal); writing -original draft (lead).Hajime Ohtsuki: Data curation (equal); formal analysis (equal); methodology (equal); writing -review and editing (equal).Shun K. Hirota: Data curation (equal); formal analysis (equal); methodology (equal); resources (equal); writing -review and editing (equal).Yoshihisa Suyama: Formal analysis (equal); methodology (equal); resources (equal); writing -original draft (equal).Jotaro Urabe: Conceptualization (lead); formal analysis (equal); funding acquisition (lead); methodology (equal); resources (equal); writing -original draft (equal); writing -review and editing (equal).
Orisarma dehaani.Ho ranged from 0.158 to 0.191, while He ranged from 0.009 to 0.258.F IS was positive in all but one population (range: 0.261 to 0.392; mean = 0.326) After removing extremely short reads and low-quality reads,17,536,469 reads (113,873 ± 2653reads per sample) were obtained from 19,215,562 raw reads (124,776 ± 2870 reads per sample).A total of 1249 SNPs across 552 loci were obtained from 154 individuals taken from 18 populations of (Table Genetic diversity parameters of the 18 O.dehaani populations in the Japanese archipelago estimated from 1249 single nucleotide polymorphisms (SNPs).
TA B L E 1Note: Significant deviation from Hardy-Weinberg equilibrium at p < .05 is denoted by asterisk.
Genetic diversity parameters of the 21 C. haematocheir populations in the Japanese archipelago estimated from 583 single nucleotide polymorphisms (SNPs).
TA B L E 2Note: Significant deviation from Hardy-Weinberg equilibrium at p < .05 is denoted by asterisk.